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It  is  suggested  here  that  in  many  envi- 
ronmental and  other  contexts  the  sever- 
ity of  an  extreme  event  might  usefully 
be  represented  by  the  sum  of  the  ex- 
cesses of  a measured  variable  over  a 
high  threshold.  The  general  form  of  the 
limiting  distributions  of  such  sums  for  a 
wide  class  of  models  has  been  derived 
by  Anderson  and  Dancy,  and  has  sug- 
gested methods  for  the  statistical  analy- 
sis of  data  concerning  extreme  severity. 
This  work  is  reviewed  here,  and  some 
extensions  to  the  distributional  theory 


are  presented.  An  application  of  the 
methods  to  atmospheric  ozone  levels, 
which  calls  for  the  extension  of  the  ap- 
proach to  take  account  of  covariate  in- 
formation is  reported. 
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1.  Introduction 

The  severity  of  a storm  or  a flood  is  often  a func- 
tion not  only  of  the  peak  value  of  whichever  envi- 
ronmental variable  is  concerned,  but  also  of  other 
aspects  of  the  extreme  event,  such  as  its  duration 
and  temporal  shape.  An  extended  run  of  days  with 
temperatures  just  below  freezing,  for  example,  can 
be  more  disruptive  to  everyday  human  activity  and 
to  animal  and  plant  life  than  a single  day  with  a 
much  sharper  frost.  Similarly,  sustained  moderately 
high  water  levels  in  a river  or  the  sea  can  lead  to 
greater  flooding  than  a more  extreme  level  lasting 
for  only  a short  time.  To  attempt  to  analyze  such 
examples  in  a way  which  captures  the  notion  of 
severity  implicit  in  them  demands  an  extension  of 
traditional  statistical  methods  for  extremes,  which 
have  tended  to  concentrate  largely  on  the  mod- 
elling of  maxima  or  storm  peaks.  In  Ref.  [1]  it  was 
suggested  that  for  an  important  class  of  applica- 
tions a simple  way  to  quantify  the  idea  of  severity 
is  in  terms  of  the  sum  of  the  excesses  of  the  envi- 
ronmental variable  over  a high  threshold  during 
the  extreme  event.  In  the  case  of  a flood,  for  exam- 


ple, this  sum  or  aggregate  excess  is  a discrete  ap- 
proximation to  the  total  volume  of  water 
overtopping  the  threshold,  and  in  the  case  of  tem- 
peratures the  analogous  quantity  defined  for  low 
values,  the  aggregate  deficit,  is  a measure  of  expo- 
sure or  cumulative  damage.  In  the  earlier  paper 
some  distribution  theory  was  developed  for  aggre- 
gate excesses,  and  an  application  to  flood  data  was 
discussed.  Here  I review  that  work  and  present 
some  extensions  of  its  distributional  results,  and 
discuss  a new  application  to  ozone  concentrations. 

2.  Preliminaries 

The  techniques  to  be  described  are  related  to 
threshold  methods  for  extremes  [2],  and  the  distri- 
butional results  are  formulated  in  terms  of  the 
Mori-Hsing  point  process  representation  [3,  4]  for 
the  structure  of  high  values  of  a stationary 
sequence.  We  briefly  recall  ideas  from  these  two 
areas. 
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Suppose  {Xj}  denotes  a sequence  of  observations, 
and  let  u be  a high  threshold.  Times  j at  which 
Xj  > u are  referred  to  as  exceedances  of  u by  {Xj}, 
and  the  sizes  of  overshoots  Xj  — u at  exceedances 
are  called  excesses  over  the  level  u.  In  environmen- 
tal applications  exceedances  are  often  found  to  oc- 
cur in  clusters  corresponding  to  physical  storms. 
Threshold  methods  are  based  on  the  modelling  of 
the  peak  excess  within  each  cluster  by  a generalized 
Pareto  distribution,  with  distribution  function  of 
the  form 

= (1) 

where  cr  > 0 is  a scale  parameter,  ^(-oo<f<<»)is 
a shape  parameter,  and  the  range  of  x is  such  that 
^xla  > — 1. 

Let  N denote  the  number  of  exceedances  within 
a cluster,  and  suppose  that  > . . . > ^/v  are  the  cor- 
responding excesses.  Then  the  suggestion  above  is 
that  the  aggregate  excess  within  a cluster 

N 

1 

is  for  some  purposes  a reasonable  measure  of  the 
severity  of  a storm  event.  For  statistical  modelling 
we  are  interested  in  the  distribution  of  S,  particu- 
larly for  high  thresholds  u.  Since  S > ^i,  we  expect 
S to  have  (in  the  limit  as  u increases)  a tail  no 
lighter  than  that  of  the  limiting  generalized  Pareto 
distribution  of  ^i.  The  distribution  of  S is  also  ex- 
pected to  reflect  the  cluster  size  and  the  pattern  of 
dependence  between  individual  excesses  . 

Suppose  now  that  M„  =maxi<;,SH  Xj.  It  is  known 
that  for  many  {X  } sequences  M„  may  be  normalized 
to  converge  in  distribution  to  some  nondegenerate 
limit.  Suppose  in  fact  that  there  is  a continuous  and 
strictly  decreasing  function  u„(t)  such  that,  for  each 

T > 0, 

lim  P{M„<u„{T))=e~\  (2) 

Let  denote  the  inverse  function  of  u„.  Consider 
now  the  two-dimensional  point  process  with  points 
(jln,u,r'(Xj)).  In  Ref.  [4],  which  generalizes  Ref. 
[3],  it  is  proved  under  a weak  long-range  mixing 
condition  A that  if  this  point  process  converges  as 
n—>oo  then  its  limit  has  points  of  the  form  (Si, 
T/Y.j),  / >1,  l<j<Ki,  where  (5,,r,),  />!  are  the 
points  of  a unit  Poisson  process  in  and  for  each 
i,  [Yjj  X}  with  y,i  = 1,  is  a point  process 


on  [1, 00 ) with  a random  number  K,  of  points.  More- 
over the  processes  [Yj-J  = I,  . . .,  A^,}  for  each  i are 
independent  of  each  other  and  of  the  {(Si,Ti)}  pro- 
cess, and  are  identically  distributed. 

A natural  interpretation  of  this  convergence  re- 
sult is  that  large  values  of  the  {Ai}  sequence  occur 
in  clusters,  located  in  time  at  the  points  of  a simple 
Poisson  process,  and  that  values  within  a cluster 
(from  the  peak  downwards)  are  given,  after  trans- 
formation, by  Ti,  TiYi2, . . . respectively  (reading  up- 
wards). Note  that,  since  the  transformation  is 
decreasing,  a cluster  peak  corresponds  to  the  lower 
endpoint  of  a vertical  string  of  points  in  the  limiting 
point  process. 

In  what  follows  it  will  be  convenient  to  suppose 
that  the  point  process  associated  with  each  cluster 
contains  infinitely  many  points  V/y  arranged  in 
increasing  order  of  size 

\ = Yn<Y.2^... 

but  that  infinite  values  of  the  Ty  are  allowed  after 
the  first  point,  so  that  K,,  the  number  of  points  in  a 
cluster,  is  just  the  index  of  the  last  finite  Yij.  By  this 
means  stochastic  properties  of  Ki  are  subsumed  no- 
tationally  in  those  of  {7/;}. 

We  are  interested  in  particular  in  clusters  of 
exceedances  by  {Xj}  of  a high  threshold  u.  Let 
v=u,r\u).  Then  Xj>u  is  equivalent  to 
w,r‘(A"y)<u,  and  so,  in  the  limit,  clusters  of  ex- 
ceedances of  u correspond  exactly  to  those  clusters 
in  the  point  process  for  which  Ti  <v.  Given  that  we 
are  dealing  with  such  a cluster  (as  we  assume  from 
now  on)  it  follows  from  the  unit  Poisson  nature  of 
{{S,,T,)}  that  T,  is  uniformly  distributed  over  (0,  v). 

For  many  {Xj}  the  transformation  u„  is  related  in 
a simple  way  to  the  marginal  distribution  function, 
F say,  of  Ay.  Suppose  in  fact  that  {Ay},  still  satisfying 
condition  A,  has  a positive  extremal  index  9.  Then 
([5],  Theorem  3.7.2) 

lim  P(M,<«„(t))=  lim  F’'Xu„ir)).  (3) 

n—*vi  /i~>* 

Hence,  if  the  tail  function  1 — F of  F is  denoted  by 
S',  it  follows  from  Eq.  (2)  that 

ndS^{u„{T)) log  F(w„(t))  ~t, 

for  large  n.  We  may  therefore  define  u„  by 

u„ir)  = S^-'(T/ne).  (4) 
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In  particular  therefore 

n9  3'{u)  = v, 

and  so  the  excesses  within  a cluster,  in  decreasing 
order  of  size,  are  in  the  limit  (dropping  the  cluster 
index  i,  no  longer  relevant) 

Cj  = u„{TYj)-u 
= 3^-^{TYjS^{u)lv)-u 
= 3^-^{T'Yj3^{v))-u, 
j = 1,  2,  ...,N  (5) 


3.1  Limit  Distributions  of  5 

Suppose  that  the  stationary  sequence  {Xj}  satisfies 
Hsing’s  mixing  condition  A and  has  positive  extremal 
index,  and  that  the  marginal  distribution  F of  the  Xj  is 
such  that  the  limiting  distribution  of  peak  excesses 
within  a cluster  is  generalized  Pareto  with  shape 
parameter  Suppose  too  that  the  corresponding  point 
process  {(jln,u,7‘(Xj))}  converges  to  a limiting  process 
with  the  structure  described  in  Sec.  2.  Then,  as  the 
threshold  level  u tends  to  the  upper  end  point,  x+  say, 
of  the  support  of  X, 

lim  p{ — 7-^  >5)  = 


where  T'  = Tlv  is  uniformly  distributed  over  (0, 1). 
The  aggregate  excess  for  the  cluster  is 

N 

s = !.{,. 

1 

where  N,  the  number  of  exceedances  in  the  cluster, 
is 


A^=max  {j:  T'Yj  <!}, 
and  T'  is  independent  of  the  Yj  process. 


3.  Asymptotic  Distributions  of  Aggregate 
Excess 

In  this  section  we  outline  various  asymptotic  dis- 
tributional properties  of  aggregate  excesses  which 
follow  from  the  preceding  discussion.  The  asymp- 
totic distribution  of  aggregate  excess  S itself  turns 
out  to  depend  on  the  Y process  partly  through  ran- 
dom sums 

«, = i z, 

/ = 1 

where  the  Z,;  are  defined  in  terms  of  {T}  by 


(^)'- 

for  ^>0 

for  ^ = 0 

-(wr 

for  ^<0. 

(6) 

for  ^>0 


( s+^\ogY\ 

""Pi 1 ) 


for  ^ = 0 


where 


for  ^ < 0 


0) 


yi{n)  = 


u 

/(i/g^(«)) 

x+  -u 


for  f > 0 
for  ^ = 0 
for  ^<0 


for  a suitable  slowly  varying  function  I,  and 


js  = min  {/■  Yj  + i^  CO  or  Rj  >s  }.  (9) 


Expectations  in  Eq.  (7)  are  taken  with  respect  to  the 
point  process  {1^}. 

This  result  is  a consolidation  and  re-statement  of 
the  main  limit  forms  found  in  Ref.  [1].  The  proof— 
essentially  a weak  convergence  argument  based  on 
the  Mori-Hsing  process  — exploits  regular  and  slow 
variation  properties  of  ^implied  by  the  assumption 
that  cluster  peak  excesses  are,  in  the  limit,  general- 
ized Pareto  distributed.  For  example,  when  ^-0, 
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^belongs  to  the  domain  of  attraction  of  the  Gum- 
bel  extreme  value  distribution,  so  that,  asx-^oo, 

~ wi{e) 

for  each  w > 0,  for  some  slowly-varying  function  / 
(see,  for  example.  Ref.  [6],  Sec.  8.13).  Thus 

c,  ~ 3^-\T’Yj3^{u))  - 


3.3  More  Explicit  Forms  for  P{Sly^(u)>s) 

When  specific  models  are  assumed  for  the  A"  pro- 
cess the  limiting  distributions  Eq.  (7)  take  on  more 
explicit  forms.  Several  examples  were  studied  in 
Ref.  [1].  Writing 

lin^  p('^_>^'j  = |(l  + sign(^)E(5,O)-‘/^for^;»^0 
/ lexp(-F(5, 0)  for^  = 0 

(11) 


(-log(r%))/(W(n)), 

asu^x+,  which  establishes  the  connection  between 
the  limiting  behaviour  of  S and  the  T-process. 

We  note  that  Eq.  (7)  reveals  in  reasonably  ex- 
plicit form  the  dependence  of  the  distribution  of  S 
on  the  number  and  pattern  of  excesses  within  a 
cluster. 

3.2  Joint  Limit  Distributions 

The  techniques  used  to  obtain  these  results  may 
be  extended  to  give  limiting  distributions  for  other 
quantities.  As  an  example  (motivated  by  a question 
from  a reservoir  engineer  about  peak  water  level 
and  total  overtopping  discharge  at  a dam  wall)  the 
joint  distribution  of  peak  and  aggregate  excesses  is 
as  follows. 

Under  the  same  assumptions  as  in  Sec.  3.1,  and 
with  the  same  notation: 


lim  P 

u—*x  + 


Vyf(M) 


>5, 


Tf  (it) 


min 

for  f > 0 

min 

for  i = 0 

min 

(“?)':  ■"-'"•I] 

for  ^<0 

(10) 


Similar  joint  limiting  distributions  may  also  be 
found  for  and  S - C\-  Like  Eq.  (10)  they  are  singu- 
lar. Methods  of  statistical  analysis  based  on  them 
have  yet  to  be  explored. 


it  was  found  that  V(s,^)  had  the  same  general  form 
in  all  cases  considered;  that  of  a concave  increasing 
function  of  s dominated  by  5 when  ^ > 0,  and  by  min 
{1,5}  when  ^<0.  See  Fig.  1. 


s 

Fig.  1.  Forms  of  the  V{s,  function. 


The  findings  and  examples  above  motivate  an  at- 
tempt to  fit  aggregate  excess  data  by  a distribution 
with  tail  function  of  this  general  form.  Two  such 
attempts  are  described  in  Sec.  4. 

3.4  Higher  Thresholds 

As  often  in  extreme  value  Statistics,  an  aim  in 
many  applications  will  be  extrapolation  to  longer 
time  periods  or  higher  levels  than  seen  in  data.  In 
particular,  for  aggregate  excesses,  extrapolations  to 
higher  thresholds  will  often  be  of  interest.  For  exam- 
ple, in  flood  applications  knowledge  of  the  aggre- 
gate excess  above  a higher  threshold  might  be  vital 
in  estimating  the  reduction  in  the  size  of  floods  that 
would  result  from  improved  river  or  sea  defences. 
The  following  presents  a simple  relationship  on 
which  extrapolation  of  aggregate  excesses  could  be 
based. 
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Suppose  that  Su  and  Su'  denote  aggregate  excesses 
above  levels  u and  u respectively,  with  u <u',  in  a 
cluster  in  which  level  u is  exceeded  (so  that  Su,  but 
not  necessarily  Su',  is  greater  than  zero).  In  a slightly 
more  refined  notation  than  used  earlier,  the  limit- 
ing forms  in  Sec.  3.1  are  limits,  'X(s)  say,  of  P{Su/ 
y^{u)>s\Su  >0)  as  u—>x+.  We  are  now  interested  in 
P{Su'lyi{u ) >s  |5u  > 0).  But 


the  general  form  of  distribution  predicted  by 
asymptotic  arguments,  and  in  particular  a simple 
Weibull  distribution  with 

P{S>s)  = t\^{-asY),  (13) 

for  some  parameters  a > 0 and  </>  was  found  to  give 
an  acceptable  fit  to  data  on  S . 


|5„  >o) 
) yd«  ) > 


P{Su'>0\Su>0) 


(12) 


for  high  u , where  is  the  peak  excess  in  the  cluster. 
Thus  the  distribution  of  aggregate  excesses  with 
respect  to  the  higher  threshold  u ' has  a point  prob- 
ability at  0 corresponding  to  the  event  /’(^i<u') 
that  no  exceedance  of  u ' occurred,  together  with  a 
form  over  the  strictly  positive  half-line  which  is  the 
the  same  as  that  of  the  original  distribution  of  ag- 
gregate excesses  except  for  an  increased  scale 
parameter.  Estimation  of  this  distribution  may 
therefore  be  based,  through  Eq.  (12)  on  estimation 
of  from  data  on  aggregate  excesses  of  u , and  of 
P{C\>u')  from  data  on  peak  excesses  of  u fitted  to 
the  Generalized  Pareto  distribution  Eq.  (1).  Rela- 
tionship Eq.  (12)  should  also  be  useful  as  a means 
of  checking  the  fit  of  specific  models  for  3^,  though 
this  aspect  has  yet  to  be  investigated. 

4.  Applications 

4.1  Floods  on  the  River  Thames 

In  Ref.  [1]  an  application  of  some  of  the  limiting 
results  above  to  data  on  levels  of  the  River  Thames 
is  described.  The  aim  was  largely  exploratory;  to  see 
whether  there  is  support  in  an  important  data  set 
for  a model  of  the  general  kind  suggested  in  Secs. 

3.1  and  3.3,  and,  if  there  is,  to  seek  an  appropriate 
parametric  form  for  the  model.  The  results  were 
surprisingly  positive:  confirmation  was  found  for 


4.2  Ozone  Concentrations 

An  analysis  of  a further  set  of  data,  which  calls 
for  the  extension  of  the  models  above  to  take  ac- 
count of  covariate  information,  is  now  reported. 

The  data  consist  of  hourly  mean  ozone  concen- 
trations at  a suburban  site  in  Stevenage,  about  25 
miles  north  of  London,  over  the  years  1978—1989. 
High  levels  of  ozone  are  known  to  cause  direct  dam- 
age to  vegetation  (see,  for  example,  Ref.  [7]).  One 
tentative  suggestion  is  that  a plant  or  tree  suffers 
damage  in  proportion  to  cumulative  exposure  to 
ozone  at  concentrations  above  some  threshold.  The 
threshold  is  not  known,  and  indeed  is  likely  to  be 
different  for  different  plants,  but  a figure  in  the 
range  40  ppb-90  ppb  might  be  plausible.  Though 
this  theory  is  at  present  no  more  than  a working  hy- 
pothesis, it  prompts  an  interest  in  the  occurrence  of 
high  values  of  aggregate  excesses  of  ozone  concen- 
trations above  moderately  high  thresholds.  The 
analysis  summarized  below  is  a preliminary  investi- 
gation into  the  possibility  of  using  the  aggregate  ex- 
cess models  of  Sec.  3 to  describe  such  high  doses.  A 
more  complete  account  of  the  biological  back- 
ground, and  of  the  application  of  the  method  to 
spatial  variation  of  exposure  over  the  UK,  is  given 
in  Ref.  [8]. 

For  the  theory  of  Secs.  2 and  3 to  be  applicable  it 
is  desirable  that  we  work  with  independent  clusters 
of  high  values.  The  hourly  data  were  therefore  sub- 
jected to  a preliminary  declustering  procedure, 
which  selected  episodes  when  concentrations  above 
a specified  ‘declustering  threshold’  were  experi- 
enced, and  ensured  that  such  episodes  were  sepa- 
rated far  enough  in  time  to  give  some  plausibility  to 
the  independence  assumption.  Figure  2 shows  a 
time  plot  of  the  resulting  aggregate  excesses  above 
a threshold  of  60  ppb,  obtained  with  a declustering 
threshold  of  50  ppb  and  with  a time  separation  be- 
tween clusters  of  at  least  48  hours  — these  values  be- 
ing chosen  as  typical  of  those  of  possible  scientific 
interest.  An  immediate  observation  from  the  plot  is 
that  the  assumption  of  stationarity  is  suspect:  the 
middle  years  1982-1986  contain  some  values  higher 
than  seen  earlier  or  later.  (There  are  known  diurnal 
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patterns  in  ozone  concentrations  too,  but  they  are 
of  too  short  a duration  to  affect  the  present  analy- 
sis.) In  view  of  the  apparent  nonstationarity  a sim- 
ple model  of  the  kind  found  useful  in  the  earlier 
analysis  would  not  on  its  own  be  expected  to  be  par- 
ticularly successful  here:  and  indeed  the  Weibull 
model  Eq.  (13)  fitted  to  aggregate  excesses  above 
60  ppb  appears  to  underestimate  the  sizes  of  the 
highest  aggregates. 

Cluster  Peaks  Aggregate  Excesses 


§ 


o - 


1978  1980  1982  1984  1986  1988  1990  1978  1980  1982  1984  1986  1988  1990 


Fig.  2.  Hourly  mean  ozone  concentrations  over  60  ppb;  Steve- 
nage 1978-1989. 

The  processes  leading  to  the  formation  of  ozone 
in  the  atmosphere  are  photo-chemical  — driven  by 
strong  sunlight.  It  is  possible  therefore  that  unusual 
weather  conditions  in  the  early  to  mid  1980s  may 
have  had  some  bearing  on  the  possible  inhomo- 
geneity. Unfortunately  sunlight  was  not  recorded  at 
the  Stevenage  monitoring  site,  nor  was  tempera- 
ture, which  is  a crude  surrogate  for  it.  Temperature 
data  were  not  readily  obtainable  either  from  nearby 
meteorological  stations,  but  were  to  hand  for 
Sheffield,  140  miles  north.  Figure  3,  showing  aggre- 
gate excesses  over  60  ppb  against  monthly  averages 
of  daily  maximum  Sheffield  temperatures,  illus- 
trates that  in  spite  of  the  geographical  separation 
there  is  nevertheless  some  connection.  It  appears 
that  the  summers  over  the  relevant  years  contained 
some  quite  warm  spells,  presumably  experienced  in 
Sheffield  as  well  as  Stevenage.  Accordingly  Weibull 
models  which  incorporate  temperature  r as  a covari- 
ate were  fitted.  Two  forms  were  used: 

P{S  >s)  = Qxp- {s/8{t)y,  (14) 

in  which  the  scale  parameter  8 depends  on  t in  the 
form  8{t)  = 5e^';  and  secondly  a model  suggested  by 


the  evidence  from  Fig.  3 that  not  all  occurrences  of 
high  temperatures  t at  the  time  of  an  ozone  cluster 
are  necessarily  associated  with  a high  aggregate 
ozone  dose.  This  suggests  a model  in  which  ozone 
clusters  are  assumed  to  be  of  two  types,  the  first 
showing  temperature  dependence  of  the  kind 
above,  and  the  second  showing  no  dependence  on 
temperature.  Thus 

f exp-(s/6(t))‘^  for  type  1 clusters 
^ 1 exp-(5/5')‘^  for  type  2 clusters 

(15) 


Monthly  average  of  daily  maximum  temperature,  degrees  Centigrade 


Fig.  3.  Aggregate  excess  ozone  over  60  ppb  vs  temperature. 

(Since  sunlight/high  temperature  is  at  best  only  one 
of  the  preconditions  known  to  be  necessary  for  the 
formation  of  ozone,  there  is  some  general  scientific 
justification  for  a model  of  this  form.)  In  fitting, 
clusters  with  aggregate  excesses  above  a specified 
level  were  taken  to  be  of  type  1.  A likelihood  ratio 
test  shows  that  model  Eq.  (15)  represents  a very 
worthwhile  improvement  over  Eq.  (14)  even  after 
allowing  for  the  inclusion  of  two  extra  parameters 
(W  = 21.26,  p < 10”\  cut-off  level  for  type  1 = 500). 
Q-Q  type  plots  for  the  two  covariate  models  are 
shown  in  Figs.  4 and  5 respectively.  (These  are  con- 
structed as  follows:  under  model  Eq.  (14)  S/{8e^') 
reduces  to  a standard  Weibull  variable  with  unit 
scale  parameter  and  shape  parameter  (f):P(SI 
{8e^')>s)  = exp(—s*).  Thus  a plot  of  the  ordered 
values  of  5/(5e^')  from  a sample  of  size  n against 
[-log(/7/j  4-1)]’^'^  should  yield  an  approximate  line 
of  unit  slope.  Figure  4 is  a plot  of  this  kind,  and 
Fig.  5 is  constructed  similarly  from  model  Eq.  (15).) 
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Fig.  4.  Q-Q  plot  for  aggregate  excess  ozone:  simple  covariate 
model  Eq.  (14). 
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Fig.  5.  Q-Q  plot  for  aggregate  excess  ozone:  two-type  covariate 
model  Eq.  (15). 

Both  plots  appear  to  show  a quite  good  fit  to  the 
Weibull  model  after  allowing  for  dependence  on 
temperature,  model  Eq.  (15)  doing  a little  better 
than  model  Eq.  (14).  Further  refinements  of  the 
models  allowing  temperature-dependence  also  of 
the  shape  parameter  <f)  gave  no  worthwhile  im- 
provement in  fit  as  judged  by  a likelihood  test. 

Though  this  is  only  a preliminary  analysis  (which 
we  hope  to  complete  with  better  temperature  data), 
the  results  so  far  are  encouraging.  They  appear  to 
show  again  that  models  of  the  form  suggested  in 
Sec.  3.3,  and  in  particular  a Weibull  model — after 
allowance  in  this  case  for  nonstationarity— can  rep- 
resent aggregate  excess  data  reasonably  well.  If  this 
is  confirmed,  then  for  example  these  models  will  be 


useful  in  estimating  return  levels  of  future  high 
doses  of  ozone  above  60  ppb  or,  following  the  re- 
sults of  Sec.  3.4,  above  higher  thresholds. 
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